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We develop one numerical method to compute the optical vortex soliton with self-defocusing Kerr- 
type nonlocal nonlinearity. With the numerical simulation method, the propagation and interaction 
properties of such optical vortex solitons are investigated. 
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I. INTRODUCTION 



Optical vortex solitons are self trapped intensity dips 
with screw phase dislocations. Owing to such screw 
phase dislocations, there is one phase singular point of 
the vortex whose real and imaginary parts both van- 
ish. Circumnavigating the phase singular point in a 
counterclockwise direction a ±2m7r phase ramp is picked 
up. The integer ±m is called a topological charge of 
the vortex. Optical vortex solitons can be generated 
in self-defocusing nonlinear medium as consequences of 
the counterbalanced effect of diffraction and nonlinear- 
ity of the medium. The numerical forms or approximate 
forms of the vortices are investigated [3, [l^ |T3|. Opti- 
cal vortices have been observed in a Kerr type nonlin- 
ear mediumfH, in a photorefractive crystal|2| and in a 
saturable nonlinear medium The stability of two 

dimensional vortices and one dimensional dark solitons 
are investigated [1 i, H, i, 0, B i, [S [ll- It is indi- 
cated that one-dimensional dark soliton stripes can de- 
cay into optical vortex solitons due to transverse mod- 
ulation instability[5, 6, J3J. A vortex of charge \m\ > 1 
is found to be topologically unstable and will split into 
m| singly charged vortices under perturbations with the 
total charge conserved 0, 0I- Hight order screw dislo- 
cation (I TO I > 1) are very sensitive to perturbations 0- 
However if the perturbation is very small, these resulted 
singly charged vortices will locate too close to be dis- 
tinguished, and in this case the splitting multicharged 
vortex is still viewed as a whole one formed as these 
nearly superposed singly charged vortices. It is indicated 
that multicharged vortices are very long-living objects 
and called to be metastable;8]. The rotation of a pair of 
first-order screw dislocation with equal signs and anni- 
hilation of a pair of dislocations of opposite signs were 
detected ji, In the limit that the interval distances 
between vortices are much larger than the size of the 
core of the vortices, the vortices can be viewed as point 
vortices [ill [l2j . Interestingly upon breakup of the input 
high-charge vortex the resulting charge-one vortices re- 
pel each other and form an array aligned perpendicular 
to the anisotropy axis of photorefractive crystal0- 



In this paper, We develop one numerical method to 
compute the optical vortex soliton with self-defocusing 
Kerr-type nonlocal nonlinearity. Such numerical method 
is similar to those presented in reference (TtI [isl . [Toj . With 
the numerical simulation method, the propagation and 
interaction properties of such optical vortex solitons are 
investigated. 



II. NUMERICAL METHOD TO COMPUTE 
NONLOCAL VORTEX SOLTIONS 

The propagation of an optical beam in spatially non- 
local self-defocusing media can be described by this fol- 
lowing (1+2) dimensional nonlocal nonlinear Schodinger 
equation(NNLSE) 



. du 

dz 
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i?(|r-r'|)|u(r',z)|2rf2r' ^ (1) 



where u(r, z) is the complex amplitude envelop of the 
light beam, |u(r, z)p is the light intensity, r = xx + j/y 
and z are transverse and longitude coordinates respec- 



tively. 



dx- 



-\- is the transverse Laplace op- 



erator. i?(|r|),(/ R{\r\)cPr — 1) is the real axisymmet- 
ric nonlocal response function, and n(r, z) — — J R{\r — 
r'\)\u{r' , z)\'^(fr' is the light-induced perturbed refrac- 
tive index. Note that not stated otherwise all integrals 
in this paper will extend over the whole transverse x-y 
plane. When i?(|r|) — 6{x,y), equation H]) will reduce 
to the local nonlinear Schodinger equation(NLSE) 



1 



(2) 



dz 2 

which has vortex soliton solutions in the form of B H, [llj 

u{r, z) = F(r)e*'^^+""^ = *(r, </J)e*^^ (3) 

where (3* = (3 is the phase constant, Lp is the azimuthal 
angle, to = ±1, ±2, • ■ • is the topology charge, F*{r) — 
F{r), F{0) — and F{r) const = 77 as r ^ oo. 

In this paper we numerically compute the vortex soli- 
ton solutions of NNLSE ([T|) expressed in the form of 
Eq. ([3]). Substituting Eq. ^ into H]) and considering 
the asymptotic behavior under r ^ cx), we have /3 
Then we have 
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77^* + ^vi* - 



i?(|r-r'|)|«'(r')|'d'r' = (4) 



FIG. 1: Modulus of vortex' amplitude F{r) with F{oo) = 
rj — l:thin solid dark line corresponds to local case w = and 
charge m = 1; thin solid red line w = 6, m = 1; dashing dark 
line w — 0,m = 2; dashing red line w = 6, m = 2; thick solid 
dark line w = 0,m — 3; thick solid red line ui = 6, m = 3. 



FIG. 2; Intensity I{x, 0, z) on plane y = of vortex of charge 
m = 2 obtained by numerical simulation with parameters 
77 = 6, w = 10. 



We discrete the function ^'(a;, y) in "ifjk = — ft- + (j — 
l)Ax, —h+{k — l)Ay^, where ~h < x < h, —h < y < h is 
the sample window, Ax — Ay is the sample step. Define 
the discrete Fourier transform (DFT) ^ by 



p,q=l 



E 

p,q=l 



(5a) 
(5b) 



where Fjk = ^exp[i2^(j - l)(k - 1)] and n 
Taking the DFT on Eq. have 
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Ax 



where Jl^fe = i ^ 

Eq. ©J we obtain 



* y" i?(|r-r'|)|*(r')pd'r' 

^)^( 
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ok 



2sin[^(fe-l)] 



(6) 



From 



{}! + ri')^,k - ^ / i?(|r - r'|)|*(r')pd2r'] 



jk 



jk 



(7) 



where /i is an arbitrary positive constant. 

We use Eq. ([7]) to iteratively compute the vortex soli- 
ton solutions. For an initial trying function, for instance, 



4'o(r, (^) = 7] 



1 — exp 



exp(im(/3). 



(8) 



where u is a constant, from Eq. ([7|), we get ^'i = .^[^"0] 
and ^'i — ^^^[^Pi]. For p > 1, we get the iteration 
scheme $p+i = ^[^p] and '^p+i = Perform 
the iteration until some accuracy is achieved, then we get 
the approximate numerical vortex soliton solutions. 

As an example, we consider this following nonlocal case 
in which the light-induced perturbed refractive index is 
governed by 



n — w'^\7^n 



(9) 



which results in 



n{r, z) 



2ttw^ 



Ko(^)|u(p,z)pdV~, (10) 



where Ko(x) is the modified Bessel function of the second 
kind and w is the characteristic nonlocal length. As indi- 
cated in Fig. ^ , the cores of vortices of the same topol- 
ogy charge and the same background intensity increase 
with the increasing of characteristic nonlocal length. The 
numerical simulation shown in Fig. Q indicates that the 
numerical solution of vortices obtained by our method 
can describe the vortices very well. As shown in Fig. ([2]) 
there is no observable splitting of the vortex of charge 
TO = 2 with absence of perturbations during the numer- 
ical simulation length z — 10. It is consistent with the 
statement that multicharged vortices are metastableQ. 

To numerically investigate the stability of singly 
charged vortex, we use an input vortex with half of the 
core size of the corresponding vortex soliton, that is 



M(r, (fi, 0) = F{2r) exp(z</?). 



(11) 



where F(r) is the modulus of the corresponding vortex 
soliton's amplitude. The simulation results are shown 
in Fig. ([3]), from which wc can find the initially shrunk 
singly charged vortex will evolve to the corresponding 
vortex soliton along with radiation of ripples. So it is nu- 
merically implies that singly charged vortex is stable. In 
another case we consider this following initial amplitude 



u{x,y,0) 



(x + d/2)-'+y^ 



xF{x,y)e''^'^'''y\ 



(12) 



which describe a field consisting of a singly charged vor- 
tex located at (x = 0, y = 0) and two normal Gaussian 
dips located at {d/2, 0) and {~d/2, 0) respectively. The 
numerical solution of Eq. ([T]) under initial condition ([12]) 
is shown in Fig. ([?]), from which, we can find that the two 
initial normal Gaussian dips without screw phase dislo- 
cation seeded will radiate ripples and become wider and 
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FIG. 3: Intensity profiles on plane y = of the evolving singly 
charged vortex with an initially shrunk core size. Dashing red 
line corresponds to input intensity of input beam with half of 
the core size of the corresponding vortex soliton; Dashing dark 
line corresponds to intensity at z = 3.3; Solid dark line z = 
8.7; Blue line z = 17.5. The solid red line is the intensity of the 
corresponding singly charged vortex soliton. The simulation 
parameters are m = 1, 77 = 6, u; = 10. 
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FIG. 4: Numerical simulation under initial condition de- 
scribed by Eq. (|12p . The left figure shows input intensity 
at = 0, a singly charged vortex located at center and two 
side-locating normal Gaussian dips initially separated by a 
distance d — 30. The right figure shows intensity at z = 7. 
The simulation parameters are w = 10, = 6, d = 30 and the 
size of each frame is equal to 27 x 27 unit. 

wider, whereas, the singly charged vortex maintains its 
shape during propagation. 

As has been previously indicated 0, that two sep- 
arated vortices of the same charge m = 1 embedded off 
axis in a finite-size Gaussian beam will rotate around the 
axis in both the linear and nonlinear cases. The rotation 
angular speed depends on the propagation distance from 
the beam waist of the background Gaussian beam and is 
not a constant. There exits a maximal rotational angle 
at infinite propagation distance from the beam waist. In 
this paper we investigate the rotation of vortices embed- 
ded in an infinite-size constant background 77. To do so, 
we consider this following initial amplitude 

uix,y,0) =-Fi{x- d/2, y)e™i'^('<-d/2,y) 

xF2{x + d/2,y)e"^^'^'-''+'^/'^-y\ (13) 

which describes two vortices of charge mi and 7712 ini- 
tially located at {d/2, 0) and {—d/2, 0) respectively. As 
shown in Figs. (O the two vortices of the same charge 
m = 1 rotate around each other. The rotational angu- 




FIG. 5: Rotation of two vortices of the same charge m — 1 
initially separated by a distance d = 4, and w = 10, rj = 6. 
The size of each frame is equal to 10 x 10 unit. 



TABLE I: The rotation period T of two vortices of the same 
charge m = 1 initially separated by a distance d. 



d 


3 


4 


5 


6 


7 


8 


T 


36.1 


50.4 


73.9 


108.8 


152.4 


200.8 


d^n/T 


0.78 


0.997 


1.06 


1.04 


1.01 


1.001 



lar speed is nearly a constant and there is no limit on 
the rotational angle of the two vortices. As indicated 
by table(|I]), the angular speed does depend on the initial 
separated distance between vortices. As has been pre- 
viously indicated [ll|, in the limit that the interval 
distances between vortices are much larger than the size 
of the core of the vortices, the vortices can be viewed 
as point vortices. For two vortices of the same charge 
TO = 1 separated by a very large distance d, based on 
the point vortices model it can be deduced that the pe- 
ripheral speed of each vortex is equal to 1/d. On the 
other hand the peripheral speed can be expressed as | ^ , 

where T is the rotation period. So we have = 1. As 
indicated by table jU, the point vortices model can give 
a rather good approximation of the peripheral speed of 
the vortices for large interval distance. 

Annihilation of two vortices of opposite charge to = 1 
and TO = — 1 initially separated by a short enough dis- 
tance are shown in Fig. It is shown that these two op- 
posite vortices annihilates each other and form a moving 
trough of the radiated ripples, which qualitatively agrees 
with the experimental observation by I. V. Basistiy, et. 
al.i7|. It is worth to note that annihilation for two far 
separated opposite vortices may be hardly possible with 
absence of other exterior interaction. Figure. ([7]) may be 
regarded as a rudimentary proof of this statement. For 
large interval distance, the point vortices model[ll|, 
predicts the directions of velocities of two vortices of op- 
posite charge to = — 1 and m — 1 are the same and per- 
pendicular to the line along these two vortices. So these 
two far separated vortices will never move close to each 
other and annihilation cannot occur. From Fig. ([7]) the 
co-moving speed of vortices is equal to 11.4/40 = 0.285, 
which is very close to 1/d = 1/4 ~ 0.25 predicted by 
point vortices model. 

Multicharged vortices are topologically unstable 0, [3]- 
In reference [2:] two possible mechanisms that split a mul- 
ticharged vortex into a set of singly charged vortices are 
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FIG. 6: Annihilation of two vortices of opposite charge m — 1 
and m — ~1 initially separated by a distance d = 3, and 
w — 10, rj = 6. The size of each frame is equal to 10 x 10 unit. 




FIG. 7: Co-moving but not annihilated two opposite vortices 
initially separated by a distance d — 4. The other parameters 
are the same as Fig. © and the size of each frame is equal to 
15 X 15 unit. 

discussed. The first mechanism is due to the fact that 
multicharged vortices are topologically unstable and sep- 
arate into a set of singly charged vortices in the presence 
of a small amount of noise, even in the framework of 
linear optics 0. The second mechanism is due to prop- 
agation effects. Anisotropic initial conditions can result 
into the splitting of multicharged vortices due to linear 
diffraction^]. Specifically the reference [2] investigated 
how an elliptically shaped high-charge vortex embedded 
in a Gaussian beam splits into charge-one vortices. In 
this paper we consider another anisotropic initial con- 
ditons that a singly charged vortex and a charge-two 
vortex initially separated by a distance co-propagate in 
a nonlocal media. As shown in Figs. ([8]) and ([9]), the 
charge-two vortex splits into two singly charged vortices 
with the present of another singly charged vortex. We 
note the nonlocality of the nonlinear response enhance 
the anisotropy of the initial conditions. In local nonlin- 
ear case, a local region's perturbed refractive index is 
solely generated by the vortex located at such a region. 
So in local nonlinear case two far separated vortices will 



never affect the perturbed refractive index of the region 
the other one located at. So the perturbed refractive in- 
dex of the region the charge-two vortex located at will 
be isotropic but not anisotropic. The charge-two vortex 
will not experience the anisotropy due to the present of 
another far separated vortex in the local nonlinear case 
and will not split into two singly charged vortices. How- 
ever in nonlocal nonlinear case, the perturbed refractive 
index of a region will depend on a distant field. The 
perturbed refractive index of the region the charge-two 
vortex located at will be anisotropic due to the present of 
another far separated vortex. Owing to such anisotropy 




FIG. 8: Interaction of two vortices of charge m — 1 and 
m = 2 initially separated by a distance of d = 4. The other 
parameters are the same as Fig. (O. 




FIG. 9: Interaction of two vortices of charge m — 1 and 
m — —2. The other parameters are the same as Fig. 



of the perturbed index, the charge-two vortex splits into 
two singly charged vortices. 

By the way, in self-focusing nonlocal case, there exist 
vortex soliton solutions which have a vanishing intensities 
rather than non-vanishing intensities in self-defocusing 
case when r — > oo. Multicharged vortices are topolog- 
ically unstable too in such self-focusing nonlocal case. 
And we will predict that the multicharged vortices will 
also split into singly charged vortics with the present of 
far separated vortices. 
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